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HUMAN TRANSFER FUNCTIONS IN MULTI -AXIS 


AND MULTI-LOOP CONTROL SYSTEMS 

By James J. Adams , Hugh P. Bergeron, 
and George J. Hurt, Jr. 
Langley Research Center 


SUMMARY 


Measurements were made of the response of a human pilot in multi -axis 
tracking tasks. These measurements are the gains in the transfer function of 
the pilot and the performance measure, root -mean- square error. The measured 
transfer functions were used to obtain analytically the closed-loop system 
characteristics. The results show that the pilot changes the response so that 
the system frequency is reduced as additional axes requiring control are added 
to the work load. It is shown that these results can be correlated with a 
theory that the pilot has a given maximum information processing capacity. 

These measured multi-axis response characteristics were used to obtain a 
quantitative description of the characteristics of a multi-loop manually con- 
trolled guidance system. The time history of the manually controlled system 
can be reproduced by using two linear analog models, one for each loop, arranged 
in series, to represent the pilot. Measurements made in the multi-axis test can 
be applied to the inner loop, and the same form of the model with modified gains 
can be used in the outer loop. It is shown that these pilot models will give an 
analytical description of the instability that can occur when there is an unex- 
pected damper failure in the system; other useful design information is also 
obtained from these models. 


INTRODUCTION 


At the present time, human response data are obtained on manually con- 
trolled closed- loop control systems through the use of simulators and the actual 
vehicle. To allow preliminary analytical design studies of flying qualities and 
to supplement pilot opinions of flying problems, an extensive effort to deter- 
mine analytical models for the human pilot is being conducted. Previous studies 
(refs. 1, 2, and 3) have measured the transfer function of human operators in 
closed-loop control systems. These studies have covered systems with a wide 
variety of controlled element dynamics in single-axis tasks, in two-axis tasks 
with one particular controlled element, and in certain nonlinear systems. The 
present study extends these measurements to include two- and three -axis tasks, 
each with two different controlled elements. As in the previous investigations. 



the transfer function of the pilot was determined by using a model-matching 
automatic -parameter tracking method. 

Another type of control situation that is encountered is the multi-loop 
problem, in which two or more variables exists each of which is dependent on the 
other. This type of situation is distinguished from the multi-axis problem, in 
which there can be two or more variables, but which are not dependent. The 
multi-loop situation is also considered in this report. The particular system 
studied represents the horizontal translation control contained in the lunar 
landing task. Time histories of the operation of such a system were obtained 
by using a fixed-base simulator, and these time histories were reproduced by 
using analog models for the pilot. Measurements obtained in the two- and three- 
axis tests were used in the inner or attitude control loop. The outer-loop 
analog pilot model was of the same form as that used in the inner loop. A 
trial-and-error method rather than the model-matching method was used to deter- 
mine the gains of this analog model. 

As an example of the use of these analog pilot models in supplementing 
pilot ratings, the multi-loop pilot transfer functions were used to study the 
unstable system characteristics that can occur when there is an unexpected 
damper failure in the system. Examples of these instabilities were obtained by 
using the fixed-base simulator, and the corresponding analytical description of 
these instabilities was also determined. 


SYMBOLS 


The units used for the physical quantities defined in this report are given 
both in the U.S. Customary Units and in the International System of Units, SI 
(ref. 4 ). The appendix presents factors relating these two systems of units. 

g gravity, 32.2 ft/sec 2 (9.81 m/sec 2 ) 

Kq_,K 2 nondimensional computer gains 

controlled element static sensitivities, either nondimensional com- 
puter gains or dimensional with units as given 

s Laplace operator, per second 

W e earth weight of vehicle 

x^ desired vehicle translation, ft (m) 

x Q actual vehicle translation, ft (m) 

X,Y,Z roll, pitch, and yaw freedom; subscript indicates total number of 
degrees of freedom 

x translation, ft (m) 
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£ damping ratio 

0 pitch attitude angle , deg 

0^ desired vehicle tilt, deg 

0 O actual vehicle tilt, deg 

t nondimensional computer gain 

cp roll attitude angle, deg 

^ yaw attitude angle, deg 

o) n undamped natural frequency, radians/sec 

A dot over a symbol denotes differentiation with respect to time, 

APPAEATUS AND TESTS 


The task performed in the multi-axis tests was compensatory tracking of a 
random disturbance. This task was performed by using a fixed-base simulator 
with a three-axis eight-ball instrument for the display, and a two-axis side- 
arm controller and rudder pedals for control. A photograph of the simulator is 

shown in figure 1. The diameter of the eight ball was 4-j- inches (11.4 cm) and 

the subject was located approximately 28 inches (71 cm) from the instrument 
panel. Fore-and-aft movement of the side-arm controller operated the pitch 
motion of the display, side-to-side controller motion operated the roll motion, 
and the rudder pedals controlled the yaw motion. The side-arm controller had a 
maximum freedom of movement of ±26° in each axis, the freedom of the rudder 
pedals was ±20°, and the sensitivity of each control was 1 volt per degree. The 
simulated controlled element dynamics used in these tests consisted of two rate 
control systems and two acceleration control systems. The transfer functions 
for the controlled element dynamics are given in the following table with maxi- 
mum accelerations and velocities (in units of instrument-display degrees) that 
will result from a 1° step displacement of the controller or rudder pedals. The 
gain of the eight ball was 5° per volt. 


Dynamics 

Maximum acceleration per 
degree stick deflection, 
deg/sec 2 

Maximum velocity per degree 
stick deflection, 
deg/sec 

0.5 

s(s + 0.5) 

2.5 

5 

2 

s(s + 1) 

10 

10 

0.5 

s 2 

2.5 

“ 

2 

s 2 

10 

! 
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These sensitivities provided adequate control for the disturbance used in 
the tests. The disturbance was obtained by passing the output of a Gaussian 
noise source through two first-order filters. The amplitude and break frequency 
for the filtered disturbance were varied during the tests to provide a suitable 
work load for each test. The break frequencies varied from 0.06 to 0.25 radian 
per second. The root -mean- square values for the disturbances are given in the 
following table: 


Break point frequency, 
rad/sec 

Root mean square, 
volts 

0.25 

9.1 

.125 

6.1 

.06 

a 4.5 


a Estimated value. 


The same disturbance signal was used on each axis of the multi-axis tests. 

The eight-ball display was used without an Euler transformation of attitude 
control torques. The task presented to the pilot was therefore one in which the 
control torque was to be applied to the gimbal axis of the instrument rather 
than to an axis system oriented to the reference marks on the instrument. At 
the same time, the disturbance characteristics were adjusted so that system 
error of less than 50 ° (except for widely spaced incidents) occurred in the 
tests. This arrangement allows the analysis of the pilot f s response to be 
treated as though each axis of the multi-axis system was a separate and 
uncoupled system. A block diagram that applies to each individual system is 
shown in sketch (a). 


Disturbance /C>v 

Di spl ay 


Pilot 


Dynamics 

<P, 0, or f 

1 

y ^ 












Sketch (a).- Block diagram of single control loop used in multi-axis simulation. 

Four NASA test pilots were used as subjects in these tests. The subjects 
were tested in the pitch, roll, and yaw axes separately, then in the combination 
of pitch and roll axes, and finally in the pitch, roll, and yaw axes combined. 

A practice run, which lasted as long as the subject desired, was given before 
each data acquisition run. Then a 3-mlmxte data rim was conducted. 

The transfer function of the pilot was obtained by the model -matching 
automatic-parameter-tracking method described in reference 1. This method uses 
a model of the form 
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Output = K 3 _t + k 1 k 2 s 
I nput ( T + s) 2 

constructed with analog computer equipment, and the gains Kp, t, and K2 are 
automatically adjusted to provide the hest possible match to the pilot. An 
illustration of the ability of this scheme to identify a known system is pre- 
sented in figure 2 . In this test the known system was of the same form as that 
selected for the adjustable model and was given gains of Kp = 7*5, t = 
and Ko = 7 - 5 - This system, or analog pilot, was placed in a control loop with 

dynamics of — - — r . The adjustable model was given initial values of gain of 

s(s + 1) 

Ki = 5, t = 10 , and Kp = 5 * The figure shows that the adjustable model quickly 
matches the known model and arrives at the desired gain values. Other examples 
in which the adjustable model is required to track time varying gains in a known 
system are presented in reference 2. 

In the multi-loop experiments, a simulator was devised to represent the 
horizontal translation guidance system, which is a part of the lunar landing 
maneuver. This system consists of a vehicle supported by a rocket alined with 
the vertical body axis of the vehicle. Translation is performed by tilting the 
vehicle until the desired horizontal acceleration is achieved. There is no 
damping of this horizontal motion in the simulation, only horizontal translation 
was studied; vertical motion and thrust control were not simulated. A pictorial 
diagram of the problem is shown in sketch (b). The simulated motion was dis- 
played to the subject with an x,y plotter and a small dial meter. The tilting 
of the vehicle was displayed on the meter, which was mounted on the moving car- 
riage of the x,y plotter. The horizontal movement of the plotter displayed 
the horizontal movement of the vehicle. A photograph of the simulator with the 
x,y plotter set upright in front of the subject is shown in figure 3. The sub- 
ject exercised control by moving the control stick from side to side. 




1000 

(305) 


Sketch (b) . - Pictorial diagram of double-loop simulation. 


K 3 

Two different inner-loop vehicle transfer functions were simulated — ^ 

Kj s 

and . The sensitivity was adjusted so that 1 degree of stick deflection 

s ( s + l) 
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produced 2.5^ deg/sec 2 of displayed tilt 

the outer-loop dynamics, which was — 

adjusted to simulate the particular task 
surface as follows: The thrust required 

the earth weight of the vehicle W 0 . It 


angle. This tilt angle is the input to 

The sensitivity of the outer loop was 

of translating a vehicle over the lunar 
to hover over the moon is one-sixth of 
was assumed that the thrust remained 


constant at this value. The expression for horizontal acceleration is then given 
hy the formulas: 


g . , 2 . , 2 

x — — sin 0 = 5-36 ft /sec sin 0 = 1.63 m/sec sin 0 

6 w e 

The small angle approximation was used to linearize this relationship, so that 
the control sensitivity of the outer loop was 

x _ 3-3^ ft/sec^ __ 1.63 m/sec 2 
0 rad rad 

Therefore the value of acceleration due to thrust for the outer loop is 

5*36 ( 1 . 63 ) when considered in terms of ft/sec 2 (m/sec^). The outer-loop simu- 
lator display sensitivity was such that 1 inch ( 2 . 5 ^ cm) of motion of the 
x,y plotter represented l 60 feet (^8.8 m) of vehicle travel. 

The problem given to the subject in the multi-loop simulation was to trans- 
late approximately 1,000 feet (305 m). The subjects operated the system so that 
this translation was completed in approximately 30 seconds. However, the tests 
were continued for approximately 2 minutes to provide a good definition of the 
hovering maneuver which followed the translation. 

The multi-loop experiments were also extended to include a third loop, which 
was displayed to the pilot by the vertical motion of the x,y plotter. Engi- 
neers only were used as subjects in the triple-loop experiments. This triple- 
loop problem is considered in the abstract sense and is not directly related to 
any physical problem. Variations in dynamics were tried in both the double - 
and triple-loop problems. Various controlled element control sensitivities were 
also investigated in the double-loop system to determine the effect of changing 
this variable. 


RESULTS AND DISCUSSION 


Multi-axis experiments .- Sample time histories taken from the multi-axis 
investigation are presented in figures 4 to 11. The measured pilot gains and 
the closed-loop system characteristics for all subjects are given in table I. 
The time histories show good agreement with the time varying results presented 
in reference 2. In the single-axis tests, the measured gains are, in general, 
constant with time whereas in the two- and three-axis tests, the gain K]_, 
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particularly in the roll and yaw axes, drops to low values or to zero for "brief 
periods of time. It is felt that these reduced values represent periods when 
the pilot has dropped all attention from that particular axis . 

Another result observed in these tests and not recognized in the results 
of reference 2 is the reduction in system natural frequency that occurs as 
additional axes for control are added to the pilot 1 s work load. These results 
are evident from the closed-loop frequencies presented in the tables. It 

can be seen that in all cases the closed-loop frequencies for a given axis are 
reduced when one or more axes are added. These changes in frequency are also 
apparent in the sample time histories. The root-mean-square values of system 
error and the normalized error, which was obtained by dividing the root-mean- 
square system error by the root -mean- square value of the disturbance, show a 
corresponding increase as the system natural frequency decreases. 

As an explanation for this reduction in system natural frequency, a theory 
is proposed that it is the effect of a maximum information processing capacity 
existing in the pilot. It is assumed that the pilot uses his maximum capacity 
in each test and when this capacity must be shared there is a reduction in the 
information processing devoted to each axis or channel. 

Shannon 1 s theorem (ref. 5 ) for channel capacity is given by the equation 
c = o> log 2 ^1 + bits per second 

where c is the information capacity of channel, co is the bandwidth of chan- 
nel, S is the average signal power, and N Is the average noise power. 

In the present investigation, the logarithmic quantity was found to be 
nearly constant. Changes in information rate showed up almost entirely as 
changes in system frequency. Since for this investigation the pilot determines 
the natural frequency of the system, it is proposed that these system fre- 
quencies correspond to the pilot's information processing capacity. 

Therefore, in the present study, the system natural frequency was con- 
sidered as a measure of channel capacity, and the ratios of system frequency 
for single-axis control to two-axis control tasks were determined for each 
subject. These results, together with the average values for all four subjects, 
are given in table II. This table shows that the system frequency for each 
variable is reduced by a factor of 0.75 when a second variable is added to the 
task. Therefore, the total capacity being exhibited by the pilot in a two-axis 
task appears to increase by a factor of 1.5 over that for a one-axis task. 

Reference 6 is a study of the information processing capacity of humans 
using one- and two-dimensional displays. This reference shows that when the 
number of dimensions of the display is doubled the ratio of information trans- 
mission capacity increases by a factor of 1.5j from 4.4 bits per second to 
6.6 bits per second. 
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The agreement between the results of the present study and that of refer- 
ence 6 suggests that a relation does exist between the maximum information 
processing capacity in a human controller and the closed-loop frequency which 
the human controller will establish in control systems of the type being 
studied in this investigation. These systems are ones in which the controlled 
element dynamics do not possess a static stability of their own. 

The ratio of approximately 0.7 between the system frequencies of single- 
and two-axis tasks also indicated that operator channel capacity can be. 
thought of as a vector which has a given magnitude and which is located in an 
orthogonal coordinate system having the same number of dimensions as the con- 
trol task. When the dimensions of the task are doubled, such as occurs when 
the task is changed from a single-axis to a two-axis task, the information 
capacity applied to each dimension will change as the orthogonal vector com- 
ponents of the original vector. Thus the capacity devoted to each channel will 

change by a factor of vl /2 or O.7O7 when going from a single-axis to a two- 
axis task. Further, when the task is changed from a single-axis task to a 
three-axis task, the capacity devoted to each dimension will change by a factor 

of v/TTF or 0.576. It can be further determined that the ratio between a two- 
axis and three -axis task should be O.815. The additional data given in table II 
Indicate that the measured data agree well with this vector concept for the 
humane information processing capacity in a closed-loop tracking task. It is 
also possible that the apparent increase in information capacity in the multi- 
axis tests resulted from the fact that the same disturbance signal was used in 
all three axes . However, the pilots upon questioning after the tests stated 
that as far as could be determined the disturbance signal appeared to be dif- 
ferent for each axis of control. 

Multi-loop experiments . - A preliminary attempt to reproduce the pilot's 
response in a multi-loop control task was presented in reference 5- This repro- 
duction was achieved by using analog models for the pilot arranged as shown in 
block diagram form in sketch (c). This scheme represents the pilot as two 
separate blocks, one In each loop, arranged In series with the loop of one 
model system closed around the other model system. An output limiter on the 
outer-loop analog pilot was also included. This limiter represents the desired 
limit of the bank angle that the pilot imposes on the control of the attitude 
of the vehicle. In reference 3,? the gains for the outer-loop analog model were 
found by a trial-and -error procedure and the gains for the inner-loop analog 
model were taken from measurements made in single-axis compensatory tracking 
tasks of reference 1. The results were sufficiently good to indicate that, in 
general, the approach was satisfactory, but it was recognized that improvements 



Sketch (c).- Block diagram of double control loop. 
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in the duplication of the manually controlled tests should he attempted. In 
particular, it was recognized that the frequency of the inner loop was too high. 
The multi-axis tests presented in the previous section indicated that, when the 
human pilot has more to attend to than just a single-axis control problem, it 
is very reasonable to consider using an analog pilot which produces a lower 
frequency to represent a human controller. 

Piloted simulator tests representative of the translation guidance portion 
of the lunar landing maneuver were made and these maneuvers were reproduced by 

using measured data taken from the multi-axis tracking studies for the inner- 
loop analog pilot. A new determination of the gains for the outer-loop analog 
pilot was also made. 

A typical time history of the piloted simulator tests and the best repro- 
duction achieved by using the analog models of the pilot are presented in fig- 
ure 12. It can be seen that the reversal in attitude angle, the small overshoot 
followed by a slow decay to the desired value of translation, and the frequency 
of the attitude mode of motion are all closely reproduced. The inner-loop 
analog pilot and associated vehicle dynamics used in obtaining the match to the 
piloted test were taken from the multi-axis test Xp of pilot D presented in 
table l(b). This particular set of measured data was chosen because the system 
frequency, which it produces, closely matched the vehicle attitude mode of 
motion frequency observed in the simulator test and because the vehicle lag 
time constant was approximately the same as that used in the simulator test. 

k 3 

The inner-loop dynamics in the simulated test was r . It is felt that 

s(s 4- 1) 

other sets of data from table I could also have been used, as Is illustrated by 
the following paragraph, or that sets of data obtained by interpolating the 
measured data could also have been used. The outer-loop analog pilot gains 
used in obtaining the match were = 0 . 09 , t = 10, and K 2 = 92 . 

Analytically determined linear system characteristics are presented in 
table Ill(a) for the values of gains used in obtaining the time history in fig- 
ure 12. Tables Ill(b) and (c) include two examples using the inner-loop pilot 
gains measured for pilot J (presented in table l(c)) and the same outer-loop 
gains used in table Ill(a). Included in table III are the closed-loop charac- 
teristics for the inner-loop alone, for the outer loop with the inner-loop 
transfer function assumed to be 1, and for the complete system. It can be seen 
from this table that some coupling does exist between the various modes of 
motion, since the characteristics of each loop considered separately are dif- 
ferent from the system considered as a whole. The damping of the attitude mode 
of motion appears to be the most sensitive characteristic because it is changed 
the most. 

The gains for the outer-loop analog model given in the previous paragraph 
provided the best match to the particular manually controlled time history pre- 
sented in figure 12. Independent variations in these gains were tried also. 

The gain Kj_ determines the frequency of the translation mode of motion. A 
50 -percent increase in K^_ resulted in an overshoot that was greater than was 
normally encountered in tests with experienced operators, and a 50 -percent 
reduction resulted in a very sluggish translation response. On the basis of an 
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overall best match to a typical manually controlled run, the value = 0.09 
was judged to be the best value. 

Values of t ranging from 2.5 to infinity were tried also. When the value 
for t was changed, K ^ and K 2 were also changed in order to keep the static 

gain K^/t and the lead time constant K 2 ^ t unchanged. A fairly good match 
was achieved with all these values of r. However, on the basis of achieving 
the best fit to both the low frequency translation mode of motion and the higher 
frequency attitude mode of motion, the value r = 10 was selected as providing 
the best overall match. 

Variations in K 2 lead to very noticeable changes in the translation mode 
of motion, as did changes, and also affected the attitude mode of motion. 

At the same time it was known that practiced human controllers would sometimes 
operate the system in a manner which resulted in a sluggish response in trans- 
lation. An example of this slower type of response is shown in figure 15 along 
with a response obtained from an analog model with an outer-loop gain K 2 of 

120. On the basis of matching the variation in responses displayed in fig- 
ures 12 and 13 , it is concluded that K 2 can vary from 9 2 to 120. 

The values of gain determined for the outer-loop analog model in this 
experiment are quite different from those measured in single-axis or multi-axis 
tracking tests. It is felt that the experiment demonstrates that a linear 
analog model of the form used, located in the system block diagram as shown, 
will adequately represent a human operator. However, in other applications, it 
is felt that the numerical values for the gains may have to be altered. 

It should be noted that the square of the frequency, or the product of the 
two lowest real roots, of the outer loop is a function of the product of the 

K 1 

static gain of the pilot transfer function and the static gain of the con- 

trolled element K^. The values for the present problem are 

2 % 

~ t“ K 3 

(0.17)(0.30) « —2(5-36) 


If the procedure presented in this report is considered for use in other 
problems that may have a different controlled element control sensitivity, then 

K-, 

it is the product Kj that should be kept constant. This adjustment is the 

first coarse adjustment that should be made in the gain to be used in the 

pilot model. Minor adjustments can then be made to suit the particular case 
under consideration. These same comments apply to the inner loop. 


Multi-loop experiments with unexpecte d changes in dynamics.- With the use 
of transfer functions for pilot response it should be possible to obtain a 
quantitative description of the instability, or pilot induced oscillation, that 
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can result when an unexpected damper failure occurs. Simulator tests were 
devised to generate this instability and the system characteristics were calcu- 
lated with the assumption that the pilot did not change his control response 
after the damper failure. Similar analyses dealing with aircraft pitch dynamics 
are presented in references 7 an( i 8* 


In the simulator tests , the previously mentioned multi-loop control task 
using the x,y plotter and dial display was used. Engineer subjects who were 

K 3 / p 

well practiced in performing this task with both — r- and K?/s dynamics 

s(s + 1) 

in the inner loop were used. They were asked to perform the task with the 

indication that the tests were to be repeat runs using — - — - dynamics in 

s(s + 1) 

k 3 

the inner loop. Then the inner -loop dynamics were changed from — to 

s(s + 1) 

Kj/s 2 during the run, simulating a damper failure. Sample results are shown 

in figure l4. The point at which the damper failure occurred is indicated in 
the figure. The typical result was the occurrence of a divergent oscillation 
or unusually large excursion in attitude angle after the damper failure. After 
5 to 10 seconds the system response returned to the type that was usually 
encountered with K^/s 2 dynamics in the inner loop. 


An analytical determination of the system instability, using the measured 
gains for pilots D and J and controlled element dynamics taken from tables l(b) 
and 1(c) for the inner-loop and the outer-loop transfer function determined in 
the previous section, was made and is presented in table III. The table com- 


pares the system characteristics with the damper in operation ^with 

2 

inner -loop dynamics for pilot D and 


0.5 


s(s + 0.5) 

. inner-loop dynamics for pilot J ) 

s(s + 1) / 

with the system characteristics with the inner-loop dynamics changed to — 


or — but with the pilot ! s transfer function unchanged. Cases using both 
s 2 

pitch and roll measured characteristics for pilot J and one case using measured 
roll characteristics for pilot D are presented. Ir can be seen that the short 
period oscillatory mode of motion changes from positive damping to negative 
damping in the tests in which the roll characteristics were used. For the case 
using pilot J pitch gains, the system did not go unstable. 


Also shown in table III are the system characteristics of the inner loop 
alone with a substitution of Kj^s 2 dynamics. These results show that the 

inner loop alone does not go unstable following the damper failure. These 
results suggest that, should an unstable condition follow a damper failure in a 
multi -loop control task, a reasonable strategy would be to cease trying to con- 
trol the outer loop and concentrate on the inner loop only. This exercise 
should continue until the pilot has become adjusted to the new dynamics. Then 
the pilot could return to controlling the multi -loop problem. The last entries 
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in tables Ill(a), (b), and (e) show the multi-loop system characteristics that 
can be expected after the pilot has become adapted to an inner -loop controlled 
element of YL^js^. These calculations show an improvement in the damping ratio 

of the attitude mode of motion, which becomes positive with pilot J but which 
remains slightly negative with pilot D. 

Effect of control sensitivity cha nges in multi -loop cont rol . - A qualitative 
investigation of the effect of control sensitivity changes in both the inner and 
outer loops on the control of the outer-loop variable was conducted. Engineer 
subjects were used. It was found that large changes could be made with little 
or no change occurring in the way in which the vehicle translation was con- 
trolled. A time history with the gains normally used in both inner and outer 
loop is shown in figure 15(a), and the effect of increasing the inner-loop 
control sensitivity by a factor of 10 is shown in figure 15(b). It can be seen 
that there is an increase in the frequency of the attitude mode of motion, but 
the translation mode of motion is unchanged. Figure 15(c) shows the result of 
leaving the inner-loop control sensitivity at normal setting and increasing the 
outer -loop sensitivity by a factor of 10. Once again the translation mode of 
motion is unchanged. The pilot still uses a limit on the attitude angle but 
this limit is reduced by a factor of 10 so that the limit value of translational 
acceleration is unchanged and, therefore, the translation maneuver is unchanged. 
Combining increases in both inner-loop and outer-loop sensitivity results in a 
small decrease in the time required to perform the translation maneuver, as is 
shown in figure 15(d). 

The combination of a lower inner-loop sensitivity and a higher outer-loop 
sensitivity also results in a negligible change in performance. A maneuver 
done under these conditions is shown in figure 15(e). The higher outer-loop 
sensitivity compensates for the low inner-loop sensitivity. As long as the 
pilot is able to control the inner loop in a linear manner (or nearly linear 
manner), as is seen in figure 15(e), there is no change in the performance of 
the maneuver. However, when the sensitivity of the inner loop was greatly 
reduced and the sensitivity of the outer loop was not increased, a condition 
resulted in which the translation maneuver was greatly slowed down. (See fig- 
ure 15(f).) The low sensitivity of the inner loop leads to a conflict between 
the pilot's desire for attitude control authority and the limit of this author- 
ity, which in the simulator was brought about by a limit in stick deflection. 
This conflict brought about a much slower control of the translation maneuver. 

Extension of the multi-loop experiments .- The multi -loop experiments, 
again using the x,y plotter and dial display, were extended to determine the 
extreme combinations of visual feedback and vehicle dyna m ics which the human 
operator could control. It was found that the subjects could control a double - 

k 3 

loop system with dynamics of — - in the inner loop. It has also been 

s^(s + 1) 

shown in references 9 and 10 that these dynamics contain the maximum controlled 
element lag with which a pilot can maintain stable control in a single-loop 
tracking task. A typical time history of the response of this double-loop 
system is shown in figure 16. The system was uncontrollable with a dynamics 

of in the inner loop. 
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It was also determined that the subjects could control a triple-loop system 
with dynamics of K-^s^- between each of the three visual feedbacks. In this 

case the three visual feedbacks were first, the dial; second, the horizontal 
displacement of the plotter; and third, the vertical displacement of the 
plotter. A block diagram of the system is shown in sketch (d). The displays 
were not meant to represent any particular vehicle motion but were presented to 
the subject as an abstract task. The subject's task was to change the vertical 
displacement of the plotter by a given distance. A typical run is shown in 
figure 17. If more controlled element lag was included in any of the three 

k 5 

loops, such as putting -5- dynamics in the inner loop, the system 

s^(s + 10 ) 

became uncontrollable. 


Second visual feedback 

First visual feedback (horizontal displacement 
(dial) of the plotter) 



Third visual feedback 
(vertical displacement 
of the plotter) 


Initial condition 


Sketch (d).- Block diagram of triple-loop simulation. 


No attempt was made to reproduce these time histories using analog models 
for the pilot. It appears that with a system as complicated as the triple-loop 
system, that linear analog models, even linear models with output limits such 
as was used in the previous sections, will not suffice. Figure 17 shows that 
the pilot used what appears to be a preprogramed type of control strategy that 
is not linearly related to any of the displayed quantities. 

It was also determined that these complex triple-loop systems were easier 
controlled with an on-off control than with the proportional control. Many of 
the records show that the pilot adapts to an on-off type of control output even 
when provided with a proportional controller. It is tentatively concluded that 
a proportional control offers no advantage in very complex and difficult to 
handle systems, and that the definite zero output position provided by on-off 
control is an advantage. The very definite dead band allows the pilot to take 
his attention away from the inner-loop display for short periods of time with 
the assurance that no inadvertent inputs will be made to the system. 

Application of pilot transf er functions to a design problem . - The pilot 
transfer functions for the multi-loop lunar landing task developed in this 
investigation were used in making certain design calculations for a full-scale 
lunar landing simulator. Certain design compromises were required in the trans- 
lation drive system of the simulator by the presence of high-frequency struct- 
ural vibrations. Sample responses to step thrust inputs for three drive con- 
figurations are shown in figure 18(a). 
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Systems 1 and 2 are simplified systems, each with a different gain setting, 
and system 3 is a more elaborate system that uses additional control signal 
inputs. The system responses define the boundary of system performance that 
can be achieved in the presence of the structural vibration problem. It was 
felt that a proper decision on the suitability of these three systems could be 
made if a calculation of system response to an input similar to that expected 
with the pilot in the loop were made. Therefore, the design analysis was 
extended to include the multi-loop horizontal translation representation 
developed in the previous section, and the response both with and without the 
simulator drive dynamics was then determined. The results are shown in fig- 
ure l8(b ) . 

The response of the vehicle alone serves as a standard with which the 
three responses, vhich include the simulator drive dynamics, can be compared. 

The results show that the differences in the three systems do affect the manu- 
ally controlled response and that system 1 gives a response that is more nearly 
like the ideal vehicle -alone response than does system 2. The apparent superi- 
ority of the alternate system, based on step input response, does not carry 
through to the manually controlled response, and therefore it was decided that 
the added mechanical complexity of the alternate system was not warranted. 
Therefore system 1 was recommended for use. 


CONCLUSIONS 


This investigation of human pilot transfer functions in multi-axis and 
multi-loop control systems supports the following conclusions: 

1. A human pilot f s response in multi -axis and double-loop control tasks 
can be well represented with linear constant -coefficient transfer functions. 

2. Measurements made in multi -axis compensatory tracking tasks show that 
the pilot changes the response characteristics, and therefore the representa- 
tive transfer function, so that the closed-loop system frequencies are reduced 
as the number of axes requiring control increases. It is shown that these 
results can be correlated with a theory that the pilot has a given maximum 
information processing capacity. 

3. Two linear transfer functions, one for each loop, arranged in series 
are used to represent a pilot ! s response in a double-loop control task. 

Measured gains from multi-axis tracking tests can be used in the inner loop. 

The gains for the outer-loop pilot model are values that correspond to a system 
frequency much lower than those found in multi-axis tests. The gains for one 
particular control task are given in this paper. 

4. These transfer functions can be used to give a quantitative description 
of system instabilities that result from unexpected vehicle dynamic changes 



such as a damper failure. Also, the pilot transfer functions obtained from the 
multi-loop simulation can he used to provide useful design information. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., October 15, 1965* 
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APPENDIX 


CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS 


The International System of Units (Si) was adopted by the Eleventh Gener al 
Conference on Weights and Measures, Paris, October i 960 , in Resolution No. 12 
(ref. 4). Conversion factors for the units used herein are as follows: 


C inches X 0.025^ = meters (m) 

Length : < 

Lfeet X 0.30^8 = meters (m) 

Acceleration: feet per second 2 X 0. 30 L 8 = meters per second 2 (m/sec^) 

Prefix to indicate multiple of unit : 10 “ 2 centi (c) 
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TABLE I.- SUMMARY OF DATA FOR MULTI-AXIS TESTS 


(a) Pilot B 


Dynamics, —7 

s(s + 1) 



Disturbance 
break point 
frequency, 
rad/sec 



Closed-loop characteristics 



Axis 

rUCttto LLI CLL &cLJ_Ilto 

Oscillatory 

Real 

roots 

Root -me an - 
square error. 

Normalized 

error 


% 

T 

*2 

rad/sec 

c 

volts 

Yl 

0.25 

37 

10 

3-5 

3.34 

0.34 

-4.75, 

- 14.0 

0.82 

0.090 

*2 

.25 

18 

7 

A 

3.31 

.36 

- 2 . 21 , 

- 10.4 

1.28 

.l 4 o 

Y 3 

.06 

7-5 

7 

5 

2.20 

.69 

- 2 . 24 , 

-9.71 

1.34 

a .298 

X 1 

•25 

33 

7-3 

2.5 

3.55 

.10 

- 3 . 27 , 

- 12.0 

1.88 

.206 

x 2 

.25 

31 

7-5 

2.5 

3.29 

.18 

- 3 - 95 , 

- 10.8 

2.4 

.263 

X 3 

.06 

10 

6.5 

5 

3.02 

.47 

- 1 . 48 , 

- 9.66 

1.34 

a .298 

Z 1 

•25 

12 

3-5 

3-5 

3.47 

.004 

- 1 . 00 , 

- 6.97 

•907 

.099 

Z 3 

. 06 

7 

3-5 

1.5 

1.94 

.08 

- 2 . 51 , 

-5.16 

1.0 

a . 222 


a From the estimated value of the disturbance root mean square. 
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TABLE I.- SUMMARY OF DATA FOR MULTI-AXIS TESTS - Continued 

Ob) Pilot D 

Dynamics. . _• disturbance break point frequency. 0.125 rad/sec 

s(s + 0.5) 


Measured gains 


K i T 


Y x 11 5.8 8.5 

y 2 15.5 5.5 4.o 
y 5 15 7.5 3.8 

X-l 25 7.5 6.5 

X 2 16 6.0 2.5 

x 5 16 7.5 3-3 

Z-L 12.5 6.0 5.8 
Z 5 5.0 10.5 5.0 


Closed-loop characteristics 
Oscillatory 


rad/sec 


Real roots 


2.70 0.21 - O . kk , -6A9 

1 . 4 l . k 2 -2.8 9, -7.41 

1.05 .37 -5.68, -9.03 

1.74 .64 -3.03, -10.23 

1.20 .26 -4.47, -7.40 

1.07 .34 -5.83, -8.94 

1.27 .62 -2.87, -8.05 

.47 .53 -9.36, -11.6 


-p. . 0.5 

Dynamics , — — 

q2 


Root -mean - 
square error , 
volts 


0.40 

.58 

1.17 

.64 

1.15 

1.60 

.84 

1.84 


Normalized 

error 


0.065 

.095 

.192 

.105 

.189 

.262 

.138 

.302 


9.5 

7.5 8.2 

0.87 0.34 

-4.97, -9.43 

O.63 

0.103 

7.0 

7.0 8.5 

.77 .34 

-4.72, -8.75 

1.54 

.252 

7.0 

8.5 8.0 

N“\ 

OJ 

I s — 

VO 

-6.62, -10.1 

1.77 

.290 

5.0 

11 10.2 

.48 .18 

-9.45, -12.4 

3.46 

.568 

3.2 

9.0 10.6 

1.58 .08 

-I.31, -4.83 

6.24 

1.02 


TABLE I.- SUMMARY OF DATA FOR MULTI -AXIS TESTS - Continued 


(c) Pilot J 


Dynamics, — disturbance break point frequency, 0.25 rad/sec 


Axis 

Measured 

gains 

Closed-loop characteristics 

Root -mean- 
square error, 
volts 

Normalized 

error 

Oscillatory 

Real 

roots 

K 1 

T 

*2 

rad/ sec 

i 

Y 1 

12 

6.0 

6.0 

3.83 

0.28 

-1.0, 

-9.83 

0.99 

0.108 

*2 

8.0 

6.5 

6.5 

3.26 

.50 

-1.0, 

- 9.76 

1.22 

.133 

Y 3 

6.0 

5.0 

5.8 

3.02 

.36 

-O.83, 

- 7.98 

l. 4 l 

.154 

X 1 

7.0 

6.5 

1^.5 

2.11 

.63 

- 2 . 25 , 

-9.06 

1.62 

b- 

H 

• 

x 2 

5.0 

5-0 

3-0 

1.72 

.47 

- 2 . 43 , 

-6.96 

2.39 

.262 

X 3 

3-0 

5-5 

5-0 

1.71 

.88 

- 1 . 5 , 

- 7.49 

2.80 

.306 

Z 1 

6.0 

5-0 

3-0 

1.93 

.42 

-2.26, 

-7.12 

1.31 

.143 

Z 3 

3.5 

3-5 

2.5 

1.71 

• 34 

-1.60, 

- 5.24 

2.27 

.249 






Dynamics, 

2 

S 2 



x l 

8.0 

7.0 

9-5 

3.32 

0.36 

- 0 . 96 , 

-10.6 

1.16 

0.127 

*2 

8.5 

7-5 

7-5 

2.52 

.47 

- 1 . 87 , 

-10.8 

1.89 

.206 

Y 3 

6.0 

8.0 

10 

2.07 

.68 

-2.0, 

-n. l 

2.5 

.274 

x 2 

5.0 

7.0 

5-0 

1-35 

.25 

- 4 . 18 , 

- 9.15 

3.41 

.383 

X 3 

5-0 

7*5 

5.0 

1.27 

.23 

-4.83, 

-9.6 

3.42 

CO 

N"\ 

ft 

Z 3 

5-5 

7.0 

7-5 

1.24 

• 35 

-5.60, 

-11.5 

3 . 7 ^ 

O 
H 
^ r 

ft 
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TABLE I.- SUMMARY OF DATA FOR MULTI -AXIS TESTS - Concluded 


(d) Pilot L 


0 5 / 

Dynamics, — r, disturbance break point frequency = 0.125 rad/sec 

s(s + 0.5) 



j — 


Closed-loop characteristics 



iyieasureu. gaina 

Oscillatory 



Root -mean- 
square error , 
volts 

Normalized 

Kl 

T 

*2 

rad/sec 

£ 

Real : 

roots 

error 

7-5 

4.0 

5.5 

1.70 

0.52 

- 0 . 89 , 

-5.84 

0.53 

0.087 

9-4 

4.8 

5-0 

1.32 

•59 

- 1 . 97 , 

-6.57 

.69 

.113 

8.0 

5.0 

5-3 

1.12 

.62 

-2.40, 

- 6.71 

.81 

.133 

7-5 

4.8 

7.8 

1.94 

.67 

- 0 . 70 , 

- 6.81 

.56 

.092 

11.5 

4.5 

2.6 

1.22 

.28 

-3.02, 

- 5.80 

1.17 

.192 

11.0 

7.5 

5.0 

1.06 

.37 

-5.43, 

-9.29 

2.5 

.41 

8.0 

3.7 

4.0 

1.43 

.38 

-1.47, 

-5.32 

.58 

.095 

6.0 

5.0 

7.4 

1.20 

.91 

-1.53, 

- 6.78 

1.12 

.183 
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TABLE III.- CLOSED-LOOP CHARACTERISTICS FOR MULTI-LOOP SYSTEMS 
(a) Pilot D, using measured pilot gains for X 2 run 
Controlled elements 


Measured gains 


Inner loop 


Ki T 


16 


Kg 


2.5 


Outer loop 


Kn 


T Kg 


Inner loop 


Outer loop 


Characteristics 


Inner loop 


0.5 


0.09 


10 


92 


s(s + 0.5) 

Outer loop 

5.56 


rad/ sec 


1.20 


Complete system 


16 


2.5 


0.09 


10 


92 


0.5 


s(s + 0.5) 


5.56 


1.11 

10.2 


i 

Real roots 


0.26 

-4.47 

-7.40 



-0.176 

-.302 

-7.60 

-11.9 


0.046 
• 99 

-0.167 

-.336 

-4.75 

-6.85 


Complete system 








16 

6 

2.5 

0.09 

10 

92 


0-5 


506 


1.20 

10.1 



Inner loop 


16 


2-5 


0.5 

S 2 


Complete system 


7-0 


8.5 


8.0 


P.09 


10 


92 


0.5 


5.36 


— • — ^ 

— . — 


1.20 

0.037 

-4.56 

-7-35 



0.737 

.197 

7.57 

10.8 

-0.094 

.98 

.99 

■ 99 
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Displayed 

error, 

volts 


Pilot 
output , 
volts 


Analog pilot 
output, volts 






0 1 2 3 

Time, min 


Figure 6.- Record of pilot J with, dynamics of 


s(s + l) 


Yaw axis of a single-axis run. 
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10 


Stick, 

volts 


Attitude angle, 
deg 


Translation, 
ft (m) 


0 


-10 


^5 

0 

■45 


0 


1000 

(505) 




Inner-loop dynamics: 


. _5_ 


s(s + 1) 


Outer-loop dynamics: 


5.36 

2 

s 



Model control 

Inner-loop model: Kj_ = 16; t = 6; Kg = 2.5; 

dynamics = ^2 

s(s + 0.5) 

Outer-loop model: - 0.09; t = 10; Kg = 120 

0 . 5.36 

dynamics = — 

s 


Figure 13.- Comparison of pilot and model in a multi-loop simulation. 


38 





(a) Control sensitivity, normal. (b) Inner-loop sensitivity, 10 x normal} 

outer-loop sensitivity, normal . 



Attitude angle, q 
deg 

A5 


Translation, 
ft |» 0 


1000 


(305) 



Time, min 


(c) Inner-loop sensitivity, normal} 
outer-loop sensitivity, 10 x normal. 


(d) Inner-loop sensitivity, 10 x normal} 
outer-loop sensitivity, 10 X no rmal . 


4o 


Figure 15*- Control sensitivity experiments. 






(e) Inner-loop sensitivity, 0.1 x normal; (f) Inner-loop sensitivity, 0.1 x normal 

outer-loop sensitivity, 10 x normal. outer-loop sensitivity, normal. 


Figure 15.- Concluded. 






Vehicle 
velocity, 
ft/sec (m/sec) 


Vehicle 
velocity, 
ft/sec (m/sec) 


Vehicle 
velocity, 
ft/sec (m/sec) 



(i) Vehicle alone , 


10'. 

(5-05) 


/ 


j 

r 

V 

/ 

\ 

/ 



r 

V 

1 

/ 

m 

: f i 

L> 

-10 ■ 

(-3.05)! 


n 

-- 









- 

- 

“1 

T 

L- 
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0 10 
Time, sec 


Vehicle 
attitude 
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(ii) Response with high gain in simulator drive mechanism. 



(iii) Response with low gain in simulator drive mechanism , 



(iv) Response with alternate simulator drive system . 


(a) 2-second step thrust input. 


(b) Calculated pilot controlled 
200-foot translation. 


Figure 18.- Comparison of response characteristics of simulator to step thrust and 

pilot controlled maneuver. 
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